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ABSTRACT 

Of the nearby transiting exoplanets that are amenable to detailed study, 
TrES-2 is both the most massive and has the largest impact parameter. We 
present z-band photometry of three transits of TrES-2. We improve upon the 
estimates of the planetary, stellar, and orbital parameters, in conjunction with 
the spectroscopic analysis of the host star by Sozzetti and co-workers. We find 
the planetary radius to be R p = 1.222 ± 0.038 i?j up and the stellar radius to be 
R± = 1.003 ± 0.027 Rq. The quoted uncertainties include the systematic error 
due to the uncertainty in the stellar mass (M* = 0.980 ±0.062 M & ). The timings 
of the transits have an accuracy of 25 s and are consistent with a uniform period, 
thus providing a baseline for future observations with the NASA Kepler satellite, 
whose field of view will include TrES-2. 

Subject headings: planetary systems — stars: individual (GSC 03549-02811) - 
techniques: photometric 

1. Introduction 

Careful follow-up observations of nearby transiting planet systems have revolutionized 
our understanding of a whole new kind of planet: hot J up iters. They have been used to reveal 
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absorption by atmospheric ato mic sodium (jCharbonneau et al.ll2002f) and the presence of an 
extended hydrogen exosphere (jVidal-Madjar et al.ll2003l ) in HD 209458b, as well as to detect 
the thermal infrared emission fr om TrES-1, HD 209458b, and HD 189733b (jCharbonneau et al 



20051 ; jDeming et al 



2005. 



2006 |). They have been used to investigate the s pin-orbit alignmen t 



of HD 209458b (IQueloz et al.ll200CH : IWinn et al.ll2005l ) and HD 189733b (IWinn et all 120061). 



Most recently, spectra of the infrared planetar y emission of HD 189733b (IGrillmair et al. 



20071 ) and HD 209458b ([Richardson et al.ll2007l ). obtained with the Spitzer Space Telescope, 
have been used to constrain models of the atmospheric content of those planets. 

Through these observations we are steadily improving our understanding of the interior 
and atmospheric structure of hot Jupiters. Future measurements, such as reflected-light 
observations or the detection of other atmospheric constituents through transmission spec- 
troscopy, will continue to advance our knowledge of these planets. One goal of the Transit 
Light Curve (TLC) Project is to support these efforts by refining the estimates of the plan- 
etary, stellar, and orbital parameters, through high-accuracy, high-cadence photometry of 
exoplanetary transits. We also seek to measure or bound any variations in the transit times 
and l ight-curve shapes that would be caused by the influence of addi tional bodies in the sys- 
tem jMiralda-EscudelhooilAgol et al.llioosl : iHolman fc Murravibood ). Along the way, we are 
exploring different techniques for photometry and parameter de termination. Previo us papers 
in this ser ies have reported r esults for t he exoplanets XO- lb (jHolman et al. 2006), OGLE- 



TR-l llb Jwinn et al.ll2007ah. TrES-1 Jwinn et all l2007bl ). OGLE-TR-lOb jHolman et al 
2007h . and HD 189733b Jwinn et ali[2O07cF " 



The present paper is concerned with TrES-2, the second tran siting hot Jupiter discovered 
by the Trans-atlantic Exoplanet Survey (jO'Donovan et al.ll2006l ). The planet orbits a nearby 
GO V star (GSC 03549-02811) and transits every ~2.5 days. Although each of the fou rteen 



known transiting exoplanets has its own story to tell, (see ICharbonneau et al.l 120071 for a 



review), the TrES-2 system has at least three distinguishing characteristics. 



(Borucki et al. 


2003; 


Basri et al. 


2005) 



Kepler will observe 

nearly six hundred transits of TrES-2 during the nominal 4 yr lifetime of the mission. This 
opportunity prompts us to improve the determinations of the orbital parameters of TrES-2 
for comparison to the future estimates from Kepler. 

Second, TrES-2 has the highest impact parameter of any known nearby transiting ex- 
trasolar planet. This makes the duration of the transit (as well as the duration of ingress 
and egress) more sensitive to changes in impact parameter. This, in turn, makes TrES-2 an 
excellent target fo r the detection of long -term changes in transit characteristics induced by 
orbital precession (jMiralda-Escudel 120021 ) . 



- 3- 



Third, the mass of TrES-2 is the largest of the known nearby transiting extrasolar 
planets. Furthermore, the radius of TrES-2 appears somewhat larger than predicted by 
simple structural models of irradiate hot Jupiters, as also appears to be the case for HAT- 
P-lb, WASP-lb, and HD 209458b (although see Burrows et al. 2007 for a contrary view). 

In what follows we present TLC results for TrES-2. In § 2 we describe the observations 
and the data reduction procedures. In § 3 we describe the model and techniques we used to 
estimate the physical and orbital parameters of the TrES-2 system, and in § 4 we summarize 
our results. 



2. The Observations and Data Reduction 



We observed four transits of TrES-2. According to the ephemeris provided by I O' Donovan et al 



(12006!) 



T C (E) = 2, 453, 957.6358 [HJD] + E x (2.47063 days), (1) 

these transits correspond to epochs 13, 15, 32, and 34 on UT 2006 Sept 11, Sept 16, Oct 28, 
and Nov 2, respectively. Observations of a fifth transit, epoch 17, were scheduled but were 
not executed due to poor weather. 

We observed these transits with KeplerCam on the 1.2m (48 inch) telescope of the Fred 
L. Whipple Observatory (FLWO) on Mt. Hopkins, Arizona. This camera (P.I. D. Latham) 



was b uilt for a photometric survey of the target field of the Kepler satellite mission (IBorucki et al. 



2003|). It has a single 4K x 4K Fairchild 486 CCD with a 23fl x 23.1 field of view. We used 
2x2 binning, for which the readout and reset time is 11.5 s and the typical read noise is 
7 e~ per binned pixel. The response of each amplifier deviates from linearity by less that 
0.5% over the range of counts from the faintest to brightest comparison star. We observed 
through the SDSS z filter, the reddest available band, in order to minimize the effect of 
color-dependent atmospheric extinction on the relative photometry, and to miminize the 
effect of limb-darkening on the transit light curve. 

The full-width at half-maximum (FWHM) of a stellar image was typically ~3 binned 
pixels (2") on Sept 11, Sept 16, and Nov 2; the FWHM ranged from ~3 to ~8 pixels on 
Oct 28. We used automatic guiding to maintain the locations of TrES-2 and its comparison 
stars to within a few pixels over the course of each night. We repeatedly took 30 s exposures 
for 3.5-5 hr bracketing the predicted transit midpoint. The conditions on UT 2006 Sept 11 
were clear during the time of the observations, and the images were taken through airmasses 
ranging from 1.05 to 1.90. The conditions on UT 2006 Sept 16 were also clear, and the airmass 
ranged from 1.05 to 1.40. There were clouds passing overhead during the observations on 



-4- 



UT 2006 Oct 28, and the airmass ranged from 1.05 to 2.50. The observing conditions were 
significantly worse during and after egress; the result was essentially observations of only a 
partial transit. Consequently, the data from Oct 28 were not included in the analysis below. 
There were very thin clouds during the observations on UT 2006 Nov 2, and the airmass 
ranged from 1.15 to 1.95. 

The images were calibrated using standard IRAfQ procedures for the overscan correc- 
tion, trimming, bias subtraction, and flat-field division. We did not attempt to correct the 
fringing that was apparent with the z filter. The fringing had a small amplitude and little 
effect on the final photometry, given the accuracy of the automatic guiding. We then per- 
formed aperture photometry of TrES-2 and 20 nearby comparison stars, using an aperture 
radius of 8.0 pixels (4'/3) for each night. We subtracted the underlying contribution from the 
sky, after estimating its brightness within an annulus ranging from 30 to 35 pixels in radius, 
centered on each star. We divided the flux of TrES-2 by the total flux of the comparison 
stars. 

To estimate the uncertainties in our photometry, we computed the quadrature sum of 
the errors due to Poisson noise of the stars (both TrES-2 and the comparison stars), Poisson 
noise of the sky background, readout nois e, a nd scintillation noise (as estimated according 



to the empirical formulas of I Young! 1 196 71 and iDravins et al.lll998l ). The dominant term is 



the Poisson noise from TrES-2. The final time series is plotted in Fig. 1 and is available in 
electronic form in Table 1. (In that table, the quoted errors have been rescaled such that 
xV^dof = 1 for the best-fitting model, as explained in the next section.) 



3. Determination of System Parameters 



Our method ology for determining the system param et ers has been described in previ 



ousTLC papers (IHolman et al.ll2006l ; IWinn et al.ll2007aU bl; iHolman et al.ll2007l ; IWinn et al. 



2007d ). and is summarized here. We assume a circular orbit of a planet (mass M p , radius 



R p ) and a star (M*, R+), with period P and inclination i relative to the sky plane. We allow 
each transit to have an independent value of T c (the transit midpoint) rather than forcing 
them to be separated by exact multiples of the orbital period. Thus, the only effect of P on 
the model is to det ermine the semimaior axis a for a given value of (M* + M p ). We fixed 



P = 2.47063 days (jO'Donovan et al.ll2006l ); the uncertainty of 0.00001 days was negligible 



1 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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Fig. 1.— Photometry of TrES-2 in the z band, using the FLWO 1.2m telescope and 
Keplercam. These data were used to estimate the planetary, stellar, and orbital parameters 
(see § 3). The bottom panel is a composite light curve created from the three data sets, 
after time-shifting and averaging into 2 min bins. The residuals (observed— calculated) are 
plotted beneath the data. 
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for this purpose. 



To calculate the relative flux as a function of the projected separation of the planet and 
the star, we employed the analytic formulas of Mandel & Agol (2002), using a quadratic 
limb darkening law, 

ii = l_ ni (l_ At )_ M2 (l_ /i )2 ? ( 2 ) 
1 i 

where / is the intensity and fi is the cosine of the angle between the line of sight and the 
normal to the stell ar surface. W e chose the values u~i = 0.22, = 0.32, based on the 



tabulated values of IClaretl (12004 ) and the estimates by ISozzetti et al.l (120071 ) of the stellar 
effective temperature, surface gravity, and metallicity. We accounted for the color- dependent 
residual airmass effects with a parameter k specific to each transit, such that the observed 
flux is equal to the intrinsic (zero airmass) flux times exp(— kz), where z is the airmass. The 
best-fitting values of k were 0.0021, 0.0086, and -0.0005, for Sept 11, Sept 16, and Nov 2, 
respectively. We also fitted for the out-of-transit flux / oot . 

The light curves cannot be used to determine both the stellar mass and radius; there is 
a fitting degeneracy R+ cx Mj . Our usual approach is to assume a value for M+ (based on 
external analyses of the stellar spectrum) and then determine R* by fitting the light curves . 
This case was slightly different because we worked in conjunction with lSozzetti et al.l (120071 ). 
who sought to improve the estimates of the stellar parameters based on the results of the 
light-curve fit. We worked iteratively, as described below in more detail; for our final analysis, 
we fixed M* = 0.98 M . 



Our fitting statistic was 



X 



N f 

£ 



/j(obs) - ij(calc) 



(3) 



where fjiphs) is the flux observed at time j, o~j controls the relative weights of the data points, 
and /j(calc) is the calculated value. It is important for o~j to include measurement errors 
and also any unmodeled systematic effects, and in particular to account for time-correlated 
noise, which effectively reduces the number of independent measurements. Our approach 
was as follows. First, we rescaled the instrumental uncertainties such that xV^dof = 1 for 
the best-fitti ng model. Table 1 lists the resulting uncertainties. Second, we followed the 
procedure of iGillon et al.l (120061 ) to decompose the observed noise into "white noise" (that 
which averages down as 1/ VN, where N is the number of data points) and "red noise" (that 
which does not average down over some specified time interval). Specifically, we calculated 
the standard deviation of the residuals (a) and the standard deviation of the time-averaged 
residuals (<Jn)- The averaging time was 1 hr (a timescale comparable to the transit event), 
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corresponding to a number N of data points that depended upon the cadence of observations. 
Then we solved for the white noise a w and red noise a r from the system of equations 



0i 



„2 



2 , 2 
<T W + <T r , 



(4) 
(5) 



Finally, to account approximately for the effective reduction in the number of independent 
data points, we rescaled the Oj in Eq. d3J) by the factor a r /(a w /y/N). In this case, the 
Sept 11 and Nov 2 transits did not show evidence for red noise according to this criterion, 
but for the Sept 16 transit the red-noise rescaling factor was 1.14. For that transit, we find 
sigma r = 0.00016 and <j w = 0.0014. To be conservative, we a pplied this same factor 1.14 
to the data from all 3 transits. 

In short there were 12 model parameters: {R*, R p , i}, as well as {T c , k, foot} for each of 3 
transits. We determined the a posteriori probability distributions for these parameters using 
the same Markov Chain Monte Carlo algorithm described in previous TLC papers. We took 
the median value of each probability distribution to be the best estimate of each parameter, 
and the standard deviation to be the 1 a uncertainty. In addition to this statistical error, 
for the special cases of and R p there is an additional error arising from the uncertainty 
in M*, which we add to the statistical error in quadrature. 

Our choice of M* merit s further discussion si nce it is based on a novel iterative procedure 
conducted in tandem with ISozzetti et al.l (120071 ). The underlying idea is that when fitting 
a light curve, the results for R+ and R p depend on the choice of M*, while the result for 

1 /3 

R+fa is independent of M* because both R+ and a vary as MJ for a fixed value of the 
orbital period. (There is, however, a minor dependence of R*/a on the choice of limb 
darkening fun ction, which is i n tur n informed by the estimates of the stellar parameters.) 



Meanwhile, as ISozzetti et al.l (120071 ) have shown, R*/a is useful for estimating M*, since it 



can be directly related t o the stellar mean density through Kepler's Third Law (see also 
Seager & Mallen-Ornelasl (bo03h ): 

1/3 

(6) 



a 

R* 



GP' X 



4vr 2 



Rl 



This makes R*/a a useful proxy for logg for the purpose of comparing the observed stellar 
properties with theoretical isochrones. The advantage of R+/a is that in typical cases it is 
more precisely determined than the spectroscopic value of log g. 



We iterated as follows. First, we fitted the light curves using the choices M* = 1.08 M , 
= 0.18, U2 = - 34, b ased on the previous estimates of the rele vant stellar paramet ers by 



O 'Donovan et al.l (120061 ). Next, we passed our results for R+/a to ISozzetti et al.l (120071 ). who 
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used it to refine the estimate of M*. (We refer the reader to ISozzetti et al.l (120071 ) for details 
on how this refinement was achieved.) In return, ISozzetti et al. (bo07h provided us with a 
new estimate of M*, along with a new estimate of the stellar surface gravity (which affects 
the choice of limb darkening law). We refitted the light curves using the updated values 
of the stellar mass and the slightly ad j usted limb darkening law. Then we passed our new 
result for R*/a back to ISozzetti et al.l (120071 ). who used it to refine the estimate of M* and 
log g, and so forth. This process converged after a few iterations, leading to the final choices 
for M*, Ui, and Ui noted above. 

While it is possible for the value of the stellar radius that minimizes x 2 to be inconsistent 
with the theoretical mass-radius relation, in this case we ha ve effectively required consistency 
with the theoretical mass-radius relation by iterating with ISozzetti et al.l (120071 ). 



Results 



The final results are given in Table [2j In addition to the results for the basic model 
parameters, we have also included in this table a number of interesting derived quantities, 
such as a/R* (which is related to the stellar mean density, as described above) and the 
calculated durations of the transit and the partial transit phases. The most interesting 
parameters are the radius of the star, the radius of the planet, the orbital inclination, and 
the mid-transit times, which we discuss in turn. 

We find the stellar radius to be R± = 1.003 ±0.027 R®, where the quoted error includes 
both the statistical error (0.017) and the systematic error due to the uncertainty in the 
stellar mass (0.021). This estimate agrees with all of th e star' s observed broadband colors 
and spectral properties as determined by ISozzetti et al.l (120071 ). as it must, given that our 
analyses were coupled as described in the previous section. 

We find the planetary radius to be R p = 1.222 ± 0.038 i?j up , where (again) the quoted 
error includes both the statistical error (0.028) and the systematic error due to the un- 
certainty in the stellar mass (0.026) . The differe n ce be tween our value, and the value 



R p = 1.220+ 

ao42 -Rjup presented by ISozzetti et al I (120071 ). is slight indeed, although our 
figure has a somewha t smal ler error bar. The reason why there is any difference at all is 
subtle. ISozzetti et al.l (120071 ) determined R p by taking our result for (R p /R±) and the asso- 
ciated uncertainty, and multiplying by their estimate for R* (which in turn was based on 
matching the observed values of T e ff, a/R*, and metallicity to theoretical isochrones). In 
contrast, we determined R p and R* simultaneously by fitting a parameterized model to the 
light curves, as described above, and then accounting for the uncertainty in the stellar mass. 
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Our analysis takes in to account the correlations between all of the parameters, while that of 
Sozzetti et al.l (120071 ) assumes (R p /R+) is independent of a/R+. In this case, our procedure 
has yielded somewhat more precise results for Rp and R±. 

For an eclipsing single-lined spectroscopic binary the surface gravity of the secondary 
(GM P / Rp, in this case) can be determined nearly independently of any assumptions regarding 
the properties of the primary(Southworth et al.ll2004l ; I Winn et al.ll2007al ; IBeattv et al.ll2007l ; 



Sozzetti et al. 2007; Southworth et al. 



20071 ). This result holds because the fitting degeneracy 



for the radial- velocity data is M p oc M^ 3 , and the fitting degeneracy for the photometric 

1 /3 

data is R p oc Mj , and in the ratio M v jR 2 v the stellar mass cancels out. There remains only 
a weak dependency of R p on the choice of limb darkening law, which is based on knowledge of 
the host star. In this case, the result is GM p /R 2 p = 1976±91 cm s" 2 , or \ogg p = 3.299±0.020. 



We confirm the finding by lO 'Donovan et al.l (120061 ) that the transit chord occurs at 
an unusually large impact parameter, b = a cos i/R* = 0.8540 ± 0.0062. This is of inter- 
est because the error in the impact parameter is much smaller when the impact parameter 
is high than when the transit is near-equatorial (all other things being equal). This fa- 
cilitates the detection of small changes in the impact parameter due to orbital precession, 
which can be caused b y additional bodies in the system or by the stellar quadrupole field 
( IMiralda-Escudel |2002| ) . A large impact parameter is also advantageous for interpreting the 
Rossiter-McLaughlin effect, as long as an acc urate external measu rement of the projected 
rotation speed of the star (vsini) is available (IGaudi &: Winnll2007l ). 



Accurate timing of exoplanetary transits is a promising method to identify additional 
planets or moons (see, e.g., Holman & Murray 2005 and Agol et al. 2005), and in this case, 
transit timing t akes on special importance because TrES-2 is in the field of view of the 
Kepler mission ( iBorucki et al.ll2003l ; iBasri et al.ll2005l ). We hav e tested whether o r not o ur 3 
measured transit times and the single transit time reported by lO'Donovan et al.l (120061 ) are 
consistent with a uniform period, by fitting a linear function of epoch number to the observed 
times. The residuals to this linear fit are shown in Fig. 2, and are consistent with zero within 
the measurement errors. Thus there is not yet any indication of timing anomalies. Based on 
our fit, we have refined the ephemeris. The new value of T c is 2, 453, 957.63479(38) [HJD] and 
the new value of the orbital period is 2.470621(17) days, where the numbers in parentheses 
are the 1 a uncertainties in the last 2 digits of each figure. 

We conducted two tests to check the robustness of our results. First, we fitted each 
of the three transits separately, and examined the scatter in the results. For each of the 
parameters {R*, R p ,i}, the 3 different results were all within the 1 a uncertainty of the 
result when fitting all the transits together. Thus the results of the 3 transits agree well 
with one another. Second, we examined the sensitivity of the results to the limb darkening 
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Fig. 2. — The timing residuals (observed — calculated) for 4 observed tra nsits, according to 
the ep hemeris of Eq. ([1]). The first point corresponds to the T c reported by lO'Donovan et al. 
( 120061 ). and the other three points correspond to the three transits reported in this paper. 
The points lie on a horizontal line, and therefore the data are consistent with a constant 
period. 
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function, finding also that the results are robust. For example, the effect on R p of allowing 
the quadratic limb- dark ening coeffici ents to be free parameters (rather than fixing them at 
the values tabulated by IClaretl (12004 )) is an increase of 1.0%. If we use a linear law instead 
of a q uadratic law, Rp is increased by 0.6%, and if we use the four-parameter "nonlinear" 
law of IClaretl (12004 ) (with coefficients fixed at ATLAS-based values) then R p is decreased 
by 0.5%. None of these changes are very significant compared to the 2.3% statistical error. 



5. Summary 



Through observations and analysis of three transits, we have improved upon the esti- 
mates of the orbital and physical parameters of TrES-2. Our resu lt s are consistent with the 
estimates of the stellar and planetary radii by lO'Donovan et al.l (120061 ). but have smaller 
uncertainties. We also show that the available transit times are consistent with a uniform 
period. In our analy sis of the photometr y we have made use of an improved estimate of 
the stellar mass from ISozzetti et al.l (120071 ). This estimate was obtained by iteratively com- 
bining values of a/R* determined from the light curves with values of effective temperature 
and metallicity determined from stellar spectra. This novel technique can be applied to all 
transiting systems for which high quality stellar spectra and high precision light curves are 
available. Our observations and analysis help lay the ground work for interpreting the ~ 600 
transits of TrES-2 that will be observed by Kepler. 
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Table 1. Relative Photometry of TrES-2 



HJD 


Relative flux 


Uncertainty 


2453989.63669 


0.9992 


0.0013 


2453989.63943 


1.0022 


0.0013 


2453989.64013 


0.9994 


0.0013 


2453989.64058 


0.9971 


0.0013 


2453989.64105 


1.0002 


0.0013 


2453989.64150 


0.9994 


0.0013 



Note. - The time stamps represent the 
Heliocentric Julian Date at the time of mid- 
exposure. The uncertainty estimates are based 
on the procedures described in § 2. We intend 
for this Table to appear in entirety in the elec- 
tronic version of the journal. A portion is shown 
here to illustrate its format. The data are also 
available in digital form from the authors upon 
request. 
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Table 2. System Parameters of TrES-2 



Parameter 


Value 


Uncertainty 


(RJR Q )(MJ0.98 M )^/3 


1.003 


0.017 


(iV-RjupXM^/0.98 Me)" 1 / 3 


1.222 


0.028 


r*/r& 


1.003 


0.027 


Rp/ Rjup 


1.222 


0.038 


R P / R* 


0.1253 


0.0010 


(W 


0.000270 


0.000012 


Mp/Mjup 1 


1.198 


0.053 


Teff [K] 1 


5850 


50 




7.63 


0.12 


z [deg] 


83.57 


0.14 


b 


0.8540 


0.0062 


irv - ti [hr] 


1.840 


0.020 


in - ti [hr] 


0.683 


0.045 


T c (13) [HJD] 


2453989.75286 


0.00029 


T c (15) [HJD] 


2453994.69393 


0.00031 


T c (34) [HJD] 


2454041.63579 


0.00030 



Adopted from Sozzetti et al. (2007). 



Note. — The system parameters and their associated uncertainties 
for TrES-2 are listed. t\, tn, and irv, correspond to the times of the 
first, second, and fourth points when the projected limb of the planet 
contacts that of the star. 



